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Abstract — Wireless networks with multiple nodes that relay 
information from a source to a destination are expected to be 
deployed in many applications. Therefore, understanding their 
design and performance under practical constraints is important. 
In this work, we propose and study three multihopping decode 
and forward (MDF) protocols for multistage half-duplex relay 
networks with no direct link between the source and destination 
nodes. In all three protocols, we assume no cooperation across 
relay nodes for encoding and decoding. Numerical evaluation in 
illustrative example networks and comparison with cheap relay 
cut-set bounds for half-duplex networks show that the proposed 
MDF protocols approach capacity in some ranges of channel 
gains. The main idea in the design of the protocols is the use 
of coding in interference networks that are created in different 
states or modes of a half-duplex network. Our results suggest 
that multistage half-duplex relaying with practical constraints on 
cooperation is comparable to point-to-point links and full-duplex 
relay networks, if there are multiple non-overlapping paths from 
source to destination and if suitable coding is employed in 
interference network states. 

I. Introduction 

One of the key technologies in next generation systems for 
achieving high throughput and providing better coverage is 
relaying. Relaying has attracted a high level of recent research 
interest with several papers focusing on various aspects of 
communicating using relays with different constraints and 
assumptions. In this work, we are concerned with the capacity 
of multistage relaying from one source to one destination 
through an arbitrary network of half duplex relays. 

An example network that we consider in detail for ease of 
explanation and clarity is the two stage relay network shown 
in Fig. [TJ In this 6-node network, the source node 5=1 
intends to communicate with the sink node D = 6 through 4 
relay nodes {Ri = 2, i?2 = 3, R3 = 4, R4 = 5} connected 
as shown. The channel gains (a, (3, 7) are shown next to the 
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Fig. 1. Two stage relay network and interference states 

corresponding edges. For simplicity, some of the gains are 
assumed to be identical. For a multistage half-duplex relay 
network such as the one in Fig. [TJ we study coding methods 



and protocols needed to achieve the best possible rate from 
source to destination for different ranges of the channel gains. 

There are two different aspects to multistage relaying when 
the relays are connected in an arbitrary fashion: (1) scheduling 
transmissions by nodes, and (2) coding methods employed 
by nodes during transmissions. One strategy for scheduling is 
to avoid interference altogether. However, the maximum data 
rate under Interference Avoidance (IA) is limited, because the 
source is nansmitting only for a fraction of the total time. 
To improve upon IA, more states of the network with the 
source node in transmit mode need to be considered. The 
scheduling task is to determine those states that are crucial 
for obtaining higher rates. When multiple nodes transmit, 
interference network states are created in the network based 
on the connectivity. Two important interference network states 
are shown in Fig. [TJ for the network of Fig. [T] In one state, 
5, R\, and R3 are transmitters and, in the other state, 5, R2, 
and i?4 are transmitters. In both the states shown in Fig. [TJ 
the source node is a transmitter and the destination node is a 
receiver. This property improves the flow of information, and 
is useful for improving the transmission rate from the source. 

In each interference network state, three different coding 
strategies of increasing complexity are considered for trans- 
mitters - Common broadcast (CB), Superposition coding (SC) 
and Dirty paper coding (DPC) for the source node alone. The 
receiving nodes in the interference network employ multiple 
access (MAC) receivers that work by successive interference 
cancellation. For different combinations of coding strategies, 
suitable rate regions are determined for each state (or interfer- 
ence network). The overall rate achievable from the source 
to the destination is computed using an optimization over 
the time-sharing of the rate regions for each state, subject to 
additional flow constraints that ensure compatibility of the rate 
vectors used for individual states. 

To place our work better, we review a sample of the 
relevant prior literature. The relay channel is a classic setting, 
introduced in (TJ, and studied extensively EJ-flU. One result 
of particular interest is the cut-set bound for half-duplex relay 
networks operating by time-sharing over a finite number of 
states [5|. This "cheap relay" bound has been used by several 
authors as an outer bound for achievable rates. 

Recently, the half-duplex diamond network with two relays 
has been studied in [6|-[9|. The multi-hopping decode and 
forward (MDF) protocol, proposed in j6|| and extended in (7), 
achieves rates close to the cheap relay cut-set bound. Wang et 
al [8 1 consider a modified diamond network with an additional 



link between the relays and propose a coding strategy using 
Dirty Paper Coding (DPC), which is shown to approach the 
cut-set bound. More protocols for general half-duplex wireless 
relay networks have been studied in ifTUl . ifTTl . 

In relation to the above, in our work, we propose and study 
multi-hopping decode and forward (MDF) protocols for a 
general relay network with half-duplex nodes in the following 
setting: (1) No cooperation is assumed for encoding and 
decoding (except in one protocol for the source node alone), 
(2) Achievable rates are compared against the cheap relay cut- 
set bound at finite SNRs, (3) The protocols and methods apply 
for a general topology of relays. The results are illustrated by 
evaluation on two different networks, where we show that the 
cut-set bound is approached for some values of channel gains. 

II. Model 

We represent a wireless network with m nodes as an 
undirected graph G = (V, E), where the vertex set V = 
{1,2, . . . ,m} represents the wireless nodes. An edge G 
E indicates that Node i and Node j are connected by an 
additive white Gaussian noise (AWGN) channel with constant 
gain denoted as h%j. 

Each node is subject to an average power constraint P and 
has a noise variance a 2 . In addition, a half-duplex constraint is 
imposed on the nodes so that they can either transmit, receive, 
or be idle at any given time. Therefore, in this work, an m- 
node half-duplex wireless network can be in M < ^# = 3™ 
states. These states are denoted Si, S2, Sm- In such a 
network, we are interested in maximizing the communication 
rate R from an arbitrary source S G V to an arbitrary sink 
D € V. The nodes in V \ {S, D} act as relays in this 
communication. Information flow from source to destination 
happens by a time-sharing of the states Sk, 1 < k < M, 
and may reach the destination in multiple hops depending 
on the connectivity of the graph. Hence, the specific problem 
considered in this work can be termed multihop, half-duplex 
relaying in an arbitrary wireless network. 

The total transmission time is normalized to one time 
unit, and state Sk is active for a fraction of the time 
(Afc could be zero) with Y^k=i ^ k = ^ ^ s m 
we assume that the state sequence and the time-sharing 
parameters are known to all nodes before transmission. Let 
Ik = {i € V : Node i is a transmitter in State Sk} be the 
set of active transmitters in State Sk, and let Jk = {i G 
V : Node i is a receiver in State Sk] be the set of active 
receivers in State Sk- When state Sk is active, simultaneous 
transmissions from nodes in Ik can interfere at one or more of 
the receivers in Jk depending on the connectivity of the nodes 
in Ik and Jk- Thus, each state Sk = (Ik, Jk) is an interference 
network lfl2l or hyperedge with Ik and Jk as the two disjoint 
vertex sets. We use the terms interference network, hyperedge 
and state interchangeably. The choice of a specific coding and 
decoding strategy for each state Sk = (Ik,Jk) determines 
possible operating rate vectors in an achievable rate region 
for that state. Since the capacity region and optimal coding 
scheme are not known for general interference networks, we 



consider three suboptimal strategies for each state based on 
different broadcast and interference processing techniques. In 
all these strategies, we impose the constraint that the receivers 
Jk cannot cooperate in decoding. Similarly, the nodes in Ik 
are assumed to encode their messages independently; however, 
in one scheme, the source is assumed to know the messages 
transmitted by the relays. 

III. Cut-set Bound 

A cut-set upper bound for half-duplex relay networks oper- 
ating by time-sharing over a finite number of states has been 
derived in 15] . This bound is presented here, briefly. 

Let and KW be the transmitted and received variables 
at node i when it is in transmit and receive states, respectively. 
The maximum achievable information rate R between source 
S and destination D in a half-duplex network is bounded as 
M 

R < sup min £ \ k I(X%y,Y$\X%), (1) 
Afc fe=i 

for some joint distributions {p(x^ 1 \x^ 2 \ ■ ■ ■ ,x^ m '\k)}, 1 < 
k < J£, where the supremum is over all Afc > such that 
Sfc=i = 1' tne minimization is over all f2 such that S G 0, 
D G n c , Xf k) = {X« : i G Cl F\ Ik}, Y$) = {Y^ : i G 
fl c n J k }, and X?" = {X^ :ieffn I k }- The above upper 
bound can be computed by solving a linear program |7*|. The 
mutual information I(X^ k y,Y^\X^) is computed exactly 
using known sum rate capacity results lfl3l when the choice 
of 51 and k results in multiple access or broadcast channels. 
When the sum rate capacity is not known exactly (e.g. for 
interference channels), the MIMO sum capacity is used as an 
upper bound. 

IV. Multihop Half-duplex Relaying Strategies 

In this section, we present the three MDF strategies that we 
propose and study in the context of a general relay network 
with half-duplex nodes. In all these strategies, the network 
operates by time-sharing between the states, where each state 
is an interference network in general. The strategies differ 
in the encoding scheme in each state. The decoder at each 
receiver employs successive interference cancellation (SIC). 

A. Common Broadcast ( CB ) Scheme 

In state Sk = (Ik,Jk), each transmitter i G Ik sends a 
common message at rate Rf? to the set of all its receivers T l _ . 
Each receiver j G Jk must decode the messages from the 
set r-j_ of all the transmitters connected to j. The decoding 
constraints at each receiver for achievability are the constraints 
for the multiple access channel corresponding to the SIC 
receiver. Therefore, the achievable rate region for each state 
Sk is defined by the following set of constraints: 

E^<^g(l+ E<e ^ P ). ( 2 ) 

for all A C and for all j G Jk- When each transmitter is 
connected to all receivers, i.e., T l _ — Jk for each i G Ik, then 
the above region is the same as the compound multiple access 
rate region in lfl4l . 



B. Superposition Coding (SC) Scheme 

In this scheme, in state Sk, each transmitter i € Ik sends 
independent messages to each of its receivers in Y l _ using 
superposition coding. Let the codeword transmitted to receiver 
j from transmitter i be x i;) -. Let the power used for this 
codeword be Pj = ctijP and be the rate. Therefore, the 
transmitter i transmits a superposition of codewords given by 
Xj = Xi 3- The received word at receiver j is 
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For simplicity of notation, we assume that the d l _ receivers 
in T l _ are arranged in descending order of channel magnitude 
from transmitter i, and Y % _ [p : q] denotes the set of elements 
of T l _ starting from the p th element to the q th element. Each 
receiver j decodes the codewords intended for itself and all 
other weaker receivers from each transmitter. Let receiver j be 
the lf L receiver in T l _ . The codewords of the weaker receivers 
r"!_[Zj + 1 : d?_\ are canceled in the SIC receiver. Therefore, 
only the codewords to the stronger receivers Y l _\\ : li — 1] 
will interfere. The received word can be written as 



Paper Coding (DPC) is used by the source to cancel inter- 
ference to its receiver caused by simultaneous transmissions 
from relay nodes. Other transmitters in Ik transmit common 
messages similar to the CB scheme. The receiver r to which 
the source is sending its DPC-coded message at rate is 
not affected by interference from other relays and will decode 
only this message. The other receivers must decode all the 
messages from all the transmitters (except the source) that are 
connected to it. For example, in the state Si shown in Fig. 
Q] S transmits a DPC-coded message to i?2 using its prior 
knowledge of the messages transmitted by Ri and R% (and 
the corresponding channel gains). Receiver R4 decodes the 
common messages transmitted by Ri and R3, and receiver 
D decodes the common message transmitted by R3. For the 
above DPC-CB scheme, the achievable rate region for state 
Sk is given by the following constraints: 
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Therefore, the achievable rate region for each state is 
defined by the following set of constraints: 



D. Optimization Model 

Now, we present the optimization problem to be solved 
to compute the achievable rate from source S to destination 
D in the multistage relay network. The optimization model 
from |fT31 is adapted to incorporate the appropriate rate region 
constraints for the MDF schemes proposed earlier. 

Let Xi A denote the information flow rate from node i to 
node j in state Sk towards the sink. Let denote the total 
information flow out of node i in state Sk- The optimization 
problem can be stated as: 



Rt < - los 
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max R, subject to: 
Flow constraints: For all i G V, we have 
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atj < 1, Vi G Ik, 
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Scheduling constraints: ^ fe Afe < 1 and > Vfc. 
Rate region constraints: The achievable rate region con- 
straints for each state depend on the encoding and de- 
coding scheme used. The rate constraints for each of the 
three proposed schemes for each state Sk are as follows: 
1) CB scheme: 



VA C Qj ^{(p,q):peT ] + ,qeT p _ [h : dt}} and Vj G J k . 

Using superposition coding allows each transmitter to send 
messages to a subset of its receivers. This receiver selection 
ability allows better spatial reuse. 

C. Dirty Paper Coding (DPC) - CB Scheme 

In the DPC-CB scheme, the source is assumed to know 
the messages transmitted by all the relays since all messages 
originate from the source. Therefore, when S G Ik, Dirty 
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where RHS of Q is the right hand side of (0. 
2) SC scheme: Equation (0|, and: 
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< A fe (RHS of ©),Vi G 4,(11) 



< A fc (RHS of ©), 



(12) 



(p,q)dA 



for all A C Qj and for all j e J^. 
3) DPC-CB scheme: 

^4 < xl Vi G 4, (13) 

^ < A fc (RHS of ©), (14) 
5^1? < A fe (RHS of ©), (15) 

for all 4cr 3 + and for all j e J k \ r. 
For the CB and DPC-CB schemes, the above optimization 
problem is a linear program. However, for the SC scheme, it 
is not a linear program since the power sharing variables a^'s 
are also optimized. Therefore, the numerical solutions for the 
SC scheme are computed using the constrained optimization 
function flnincon in MATLAB. 

V. Numerical Results 

We evaluate and compare the rate achieved by the MDF 
schemes: (1) CB, (2) SC, and (3) DPC-CB for two different 
network topologies and channel realizations. The cheap relay 
cut-set upper bound for half-duplex relay networks and the 
rate achieved by the IA scheme are also evaluated. The rate 
achieved by each scheme is obtained by solving the optimiza- 
tion problem in (0 with appropriate rate region constraints. 

Since the diamond network has been studied in detail in Q, 
0, we skip details and simply mention that the proposed MDF 
protocols recover similar results for the diamond network. 

A. Two stage relay network 
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Fig. 2. a = 1, 7 = 1, vary 0. 

We first consider the two-stage relay network shown in Fig. 
Q] For evaluating the cut-set bound, all the 2 2 • 3 4 = 324 states 
were considered. The states that avoid interference (called IA 
states) are the states with a single transmitting node. For the 
proposed MDF protocols, interference network states with two 
transmitters (QM = 10 states) and some states with three 
transmitters (5 out of ( 3 ) = 10 states) are used along with 
the IA states. Two of the states with three transmitters are 
shown in Fig. Q] for illustration. In Fig. [2] the cut-set bound, 
determined by the source cut, is at 1 for all f3. For large 



j3, the states used are S3 = ({S, R2, R3}, {Ri, R4, D}) and 
S4 = ({S, Ri,Ri}, {i?2, R3, D}). The receivers in both these 
states see strong interference, which can be canceled at the 
receiver. For instance, in state S3, the receiver Ri can decode 
the source's message in the presence of strong interference 
from i?2 and R3. Because of this, all three MDF schemes 
achieve capacity of 1 by equal time-sharing of states S3 and 
S4. For small /?, common broadcast at the relays is limited 
by a weak receiver with close-to-zero capacity. Superposition 
coding, which enables different rates to receivers, proves to be 
better at low values of (3. For SC, states Si and S2 (shown in 
fig- [H> are chosen, and the rate is limited by the interference at 
relays i?i and i?2- DPC is marginally weaker, since the relays 
continue to do common broadcast when the source does DPC. 
However, when /3 = 1 (0 dB), DPC is better as SC becomes 
identical to CB for identical channel gains. 
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Fig. 3. a = 1, /3 = 1, vary 7. 

In Fig. [3] the cut-set bound, determined by the source cut, 
is at 1 for all 7. For 7 > 1, the DPC-CB scheme achieves 
the cut-set bound for lower 7 than SC and CB schemes. DPC 
achieves capacity by time-sharing the states Si and S2. The 
interference at relays Ri and R2 are canceled using DPC, 
while the interference at R3 and R4 is overcome because the 
gains of the R2 — > R3 and R\ — > R4 links increase with 7. 
The same states are used for the SC scheme as well. However, 
interference at Ri and i?2 are overcome only for larger 7. For 
very large 7, CB scheme also approaches the cut-set bound 
by time-sharing between the states ({S, R4}, {Ri, D}) and 
({i?i}, {-fiU})- For small 7, we see that DPC-CB achieves a 
rate of 0.7, while the SC and CB achieve rates of 0.67 and 
0.55 respectively. For both DPC-CB and CB schemes, states 
S3 and S4 are chosen. While the interference at R3 and R4 
limits the DPC-CB scheme, the CB scheme is limited by the 
interference at relays R\ and i?2- 

B. Rectangular grid network 

Consider the 4x3 rectangular grid network shown in Fig. [4] 
Since the number of possible states is prohibitively large, we 
first select three non-overlapping paths from the source node 



S = 2 to the destination node D = 11. We know from the 
two-stage relay network example that multiple flow paths used 
appropriately with interference processing can be effective. 
The paths chosen are # ^ 4 -> 7 L>, S" 5 8 D 
and S — > 6 — 5- 9 — > D. Using the nodes on these paths, the 
three states chosen for scheduling are ({S, 6, 8}, {4, 9, D}), 
({S, 4, 9}, {5, 7, D}), and ({S, 5, 7}, {6, 8, D}). Note that the 
source node is a transmitter and the destination node is a 
receiver in all three chosen states. Also, the other two trans- 
mitters are chosen to be at different distances from the source. 
With this choice of states, we have a two-stage relay network 
with six relay nodes {4,5,6,7,8,9} aiding communications 
from the source to the destination. 




Fig. 4. 4x3 Grid Network. 

In Fig. [5] the gains f3 and 7 are set to 1, and the gain a 
is varied. We notice that the DPC-CB scheme approaches the 
capacity for a large range of values of a > 1. The CB and 
SC schemes are limited by the interference at relays 4 and 5 
even for large a. For small a, the DPC-CB and CB schemes 
are limited by the common broadcast constraint at the relays. 
While SC scheme can perform better, it is still limited by 
interference at relays 4 and 5 compared to the cut-set bound. 
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Fig. 5. Performance in Grid Network, = 1,7= 1, vary a. 

In summary, in larger networks, the choice of schedule is 
important. We have used a path-based heuristic and relied on 
interference-processing for approaching the cut-set bound. 

VI. Conclusions 
Based on this work, two interesting comparisons are pos- 
sible for multistage half-duplex relay networks based on the 



cut-set bound. For the network in Fig. [T] the cut-set bound 
evaluates to C vv = log(l + a 2 P/a 2 ), which can be interpreted 
as the capacity of a point-to-point link with power constraint 
P and channel gain a. Using the protocols in this work, we 
have shown that rates up to C vv are achievable by multistage 
half -duplex relaying in the network of Fig. [T]for certain ranges 
of the channel gains a, /3 and 7. A necessary condition for 
achieving the point-to-point capacity under the half-duplex 
constraint is that the source needs to be in transmit mode at all 
times. From our work, it appears that continuous transmission 
by the source and information transfer through the half-duplex 
relays is possible as long as there are two or more non- 
overlapping paths from the source to the destination (which is 
true in Figs. Q] and 2]). Further, coding in interference networks 
created by multiple transmitters and receivers of the relay 
network is crucial for enabling the information flow. 

The second comparison is with full-duplex relays. The 
achievable rate even with full duplex relays is bounded by the 
sum rate across the source-broadcast cut, which is equal to 
C pp , for the network in Fig. Q] Once again, we observe from 
our work that two non-overlapping paths through the relays 
and interference-network coding enable a half-duplex relay 
network to achieve the full-duplex cut-set bound for certain 
ranges of channel gains. 
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